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Several recent studies have dealt with the effects of short range correlations on the momentum 
distribution of nucleons in nuclei. Here we investigate the correlation effects on the density and 
spectral distribution in coordinate space. A combination of experimental data and spectral functions 
calculated from realistic N-N interactions allows us to resolve a recently uncovered discrepancy with 
occupation of quasi-particle states derived from (e,e'p) data. 
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I. INTRODUCTION 

Much of the understanding of atomic nuclei is based on 
the assumption that nucleons move, independently from 
each other, in the average potential created by the inter- 
action with all other nucleons. A more fundamental ap- 
proach to the understanding of nuclei has to start from 
the underlying nuclcon-nucleon (N-N) interaction. Re- 
alistic models of the N-N interaction exhibit a strongly 
repulsive central interaction at small inter-nucleon dis- 
tances and a strong tensor component. These features 
lead to properties of nuclear wave functions that are be- 
yond what is describable in terms of Independent Particle 
(IP) motion. In particular, strong short-range correla- 
tions (SRC) are expected to occur. 

The effects of the short-range correlations are known 
for systems where an accurate solution of the Schrodingcr 
equation for a realistic N-N interaction can be obtained 
Very light nuclei (today up to A<10) and infinite 
nuclear matter are amongst the systems where this is 
feasible 0, 0, . The corresponding calculations show 
that in a microscopic description of nuclear systems the 
short-range and tensor parts of the N-N interaction have 
a very important, not to say dominating, influence with- 
out which not even nuclear binding can be explained. 

Due to these short-range correlations the momentum 
distributions of nucleons acquire a tail extending to very 
high momenta k and part of the strength, located in IP 
descriptions at low excitation energy E, is moved to very 
high excitation energies. 

In the past, most experimental investigations were con- 
fined to rather low momenta and energies, i.e. to the 
region where the strength is dominated (but not entirely 
given) by the IP properties. In this region, the conse- 
quences of short-range correlations are indicated primar- 
ily by a depopulation of states in comparison to the pre- 
dictions of IP models (including the long-range correla- 
tions which can be described by configuration mixing). 
According to the calculations mentioned above, a depop- 
ulation of the order of 20% is expected. 

/From the experimental information available up to 



now, the depopulation of IP strength at low k, E is es- 
tablished (for a caveat see below) . Much less is known 
from a direct measurement of the strength of the spec- 
tral function S(k, E) at large k and E. A recent (e,e'p) 
experiment, performed at high momentum transfer q in 
parallel kinematics by Rohe et al., provides the first di- 
rect measurement |^|^. 

This correlated strength has always been discussed in 
k, -E-space where (part of) it can be separated from the 
IP strength. In this Rapid Communication, we take an 
orthogonal look at the correlated strength and discuss it 
in coordinate space (r-space). We address this question 
from both the theory and experiment side. 

The trigger for this study lies in difficulties experienced 
in the past in interpreting data in terms of IP models. For 
example, fits with IP wavefunctions of the nuclear charge 
density often yielded form factors (of IP dominated tran- 
sitions) with incorrect g-dependence. Fits with IP wave 
functions also have difficulties to reproduce the total den- 
sities in the nuclear interior. The origin of these difficul- 
ties: total densities have contributions from correlated 
nucleons that do not appear in observables dominated 
by quasi-particle properties. The correlated nucleons pre- 
sumably have a different radial distribution. 

The goal of this paper is to derive quasi-particle (QP) 
and correlated distributions in r-space using Green's- 
function theory. We compare the results to the correlated 
density in r-space which we obtain from the difference of 
the density — known from elastic electron scattering — 
and the QP contributions known from (e,e'p) reactions. 

As a side- product, this study also sheds light on a re- 
cently uncovered problem with QP occupation numbers 
derived from (e,e'p) experiments with low and high q, 
respectively Q. 



II. SINGLE-PARTICLE SPECTRAL FUNCTION 

The evaluation of the single-particle spectral functions 
for 12 C, the nucleus we use for our study, has been per- 
formed within the framework of the Green's function 
method Q using the techniques described in 0, ITot Hll| . 
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The nucleon self-energy Ey (p m ,p n > E) is determined in a 
discrete basis of Bessel functions 4> Pm ij(r) with appropri- 
ate boundary conditions at the surface of a spherical box 
with radius R\) 0X - These basis states are identified by 
the angular momentum quantum numbers I and j and a 
radial quantum number (p m , Pn) referring to the momen- 
tum. For a box radius R\^ ox of typically 20 fm, it turns 
out to be sufficient to include around 60 basis states for 
each partial wave. 

The self-energy contains a Hartree-Fock contribution 
Ey F calculated in terms of a nuclear matter G-matrix 
plus complex correction terms, A£y , which account for 
the inclusion of two-particle one-hole and two-hole one- 
particle contributions. These correction terms are calcu- 
lated directly for the finite nucleus 12 C, describing the in- 
termediate particle states by plane waves orthogonalizcd 
with respect to the occupied hole states. This is a good 
approximation to describe the effects of SRC, however, 
it tends to underestimate the spectral strength due to 
long-range correlations at missing energies slightly above 
the two-hole one-particle threshold. 

The single-particle Green's function is determined from 
this self-energy + A£y by solving the Dyson equa- 
tion in the box basis described above. From the imagi- 
nary part of this Green's function one can calculate the 
spectral function in this basis |lll | or determine it in con- 
figuration space by the transformation 

Sij{r,r';E) = ^ Pm y (r)Sij (p m , p„ ; E)<fi* n l j (r' ) , 

using the Bessel functions 4> Pm ij(r) described above. The 
spectral function can be split into the QP contribution 
Sy , which only occurs in the s 1 / 2 and p 3 / 2 partial 
waves, and in the continuum contribution, Sf° nt , which 
originates from the imaginary components in the self- 
energy. This leads to the single-particle density 

dESf° nt (r,r-E) 

= PQp(r) + Pcorr(r) , 

where the integration over missing energies E starts at 
the threshold of two-hole one-particle configurations. We 
have assigned the label corr to the part of the single- 
particle density, which originates from the continuum 
part of the spectral function to indicate that this cor- 
related density is absent in the IP model. 

Results for these contributions to the point density of 
protons in 12 C are displayed in Fig.QJ The calculation of 
these densities have been performed using the CD-Bonn 
potential for the N-N interaction ^3] • A fraction of the 
proton density, which accounts for around 5 protons, is 
described by the QP part and the rest is covered by the 
correlated density. 

In order to allow for a better comparison of the radial 
shape of the density contributions, Fig. ^also contains 
the correlated density p corr (r) multiplied by a factor of 3. 



The comparison shows quite clearly that the correlated 
density is located much more in the center of the nucleus 
than the QP contribution. 

The correlated single-particle density is distributed 
over partial waves Ij including those which are unoc- 
cupied in the IP model. A large fraction of the cor- 
related strength, however, is contained in partial waves 
with I = (around 31%) and I = 1 (around 37%). One 
also should note that, contrary to what one naively could 
expect, the strength in the higher I states does not con- 
tribute at large r; the corresponding large values of E 
pull the radial wave functions to lower r. 
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FIG. 1: Comparison of QP and correlated densities from the- 
ory. 



III. DENSITIES FROM (e, e) 

For 12 C an extensive set of elastic electron scattering 
data is available [l3l ITiL ITEL ITfij covering the range of 
momentum transfer between 0.13 and 3.7 fm . The 
Carbon rms-radius is precisely known from a /x-X-ray 
experiment |T^ . 

These data have been used to determine a model- 
independent charge density using the SOG-approach of 
[T^j . The procedure employed has been described in 0] . 
This yields the charge density as a function of r, together 
with an error bar that covers both the statistical and sys- 
tematic uncertainties of the data, as well as the model- 
error. 

In order to obtain the point density, we have unfolded 
the effect of the finite size of the proton and neutron 
charge density. This has been done by parameterizing 
the point density, folding it with p p (r)+p n (r) and fitting 
the resulting density to the charge density as determined 
above. The contribution of the electro-magnetic spin- 
orbit term to the charge density turned out to be neg- 
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ligiblc. The resulting folded density agrees within ~1% 
with the one given in [20l |. 

The resulting point density is shown in figs. 12131 The 
error bars in general are too small to be seen. 

IV. QP ORBITS FROM (e, e'p) 

For Carbo n, q uite an extensive set of (e,e'p) data is 
available [U H£ H |H |p ; a compilation 

is discussed in Part of this data has been taken at 
low q with the goal to determine the lp and Is quasi- 
particle momentum distributions and occupation num- 
bers. Some data have been taken at large q mainly in 
connection with the determination of nuclear transparen- 
cies for high-energy protons. 

The low-q data, taken with good energy- and mo- 
mentum resolution mainly at NIKHEF and Saclay, have 
been analyzed in the framework of DWBA using optical 
potentials known from proton-Carbon scattering. The 
QP radial wave functions have been parameterized using 
Woods-Saxon (WS) potentials. Lapikas et al. Q have 
made a coherent analysis of the entire data set. The 
occupation of the QP orbits, obtained by summing the 
experimental spectroscopic factors, turns out to be rather 
low in comparison to what is known for other nuclei Q; 
the summed lp plus Is strength amounts to 3.4 protons 
only (56% occupation). 

The high-g data, determined in part with moderate 
energy- and momentum resolution, were taken mainly at 
SLAC and JLAB. The data were summed over a large 
region of initial momentum k and removal energy E, 
and fitted using WS radial wave functions and theoreti- 
cal transparencies. When using the most reliable trans- 
parencies, confirmed by other experiments, Lapikas et al. 
found a much higher occupancy of the QP orbits, 5.0-5.6 
protons (^87% occupation). 

This discrepancy — which is very embarrassing to the 
practitioners of (e,e'p) as it sheds serious doubts on the 
quantitative interpretation of (e,e'p) data — obviously 
needs to be better understood, and is discussed in more 
detail below. This difference also has led to speculations 
about g-dependent QP-occupations, for which we see no 
physical basis. 

A partial reason for the difference between the low- 
q and high-g results is immediately clear: The low-g 
data cover the region of missing momenta of typically 
<180MeV/c and missing energy E <50MeV, the high-g 
data extend to 300MeV/c. The high-g data also cover a 
larger range in missing energy, they are integrated up to 
typically 80MeV. In this larger k, E-i&nge, there is not 
only QP strength, but also a fraction of the correlated 
strength is integrated over. Before making a valid com- 
parison, this correlated strength needs to be removed. 

In order to correct for this effect, we start from the 
high-g (e,e'p) data taken in a recent JLAB experiment 
in quasi-elastic kinematics, which minimize final state 
interactions FSI and meson exchange current contribu- 



tions MEC. This experiment yields, in agreement with 
the previous JLAB and SLAC experiments, 5.2 protons 
in the integration region E<80MeV, k<300MeV/c. This 
number we correct for the continuum contribution using 
the calculated spectral function discussed above. With 
this correction the discrepancy between the low-g and 
high-g results is significantly reduced; the QP occupa- 
tions now are 3.4 vs 4.5 protons for the low-g/high-g data, 
respectively Given the uncertainties of these numbers — 
believed to perhaps 10% — there still is a worrisome in- 
compatibility. 

In order to proceed, we have to choose. We have de- 
cided to use the occupation coming from the high-g mea- 
surements, as we judge the interpretation of these data to 
be more safe. The low-g data suffer from uncertainties 
in the treatment of the final state interaction. Due to 
the low energies of the outgoing proton (70MeV for the 
NIKHEF data [25^ coupled channel effects not treated 
in the usual DWBA analysis should be relevant. Van der 
Steenhoven et al [25l l29j have shown that inclusion of 
these effects would increase, for the rather soft nucleus 
12 C, the QP occupation by ^20%. For the kinematics 
of the low-g experiments, the calculations of Boffi et al. 
[30l l3l| also predict significant MEC effects that would 
lead to a further increase of the QP occupation. 

The value for the QP occupation deduced from the 
high-g measurement is also compatible with the corre- 
lated strength measured directly in the recent experiment 
by Rohe et al. [fj. This measurement agrees with the- 
oretical predictions for the correlated strength of ~20%. 
The summed QP strength from the low-g (e,e'p) data 
(which thus includes the fragmentation due to long-range 
correlations), on the other hand, would correspond to 
>40% correlated strength, i.e. be unrealistically high. 
Furthermore, as we will see below, our choice of the QP 
occupation is confirmed by a consistency check in our 
analysis. 

Before proceeding to the calculation of the QP density, 
one more effect of correlations must be removed from the 
results of Lapikas et al. While the p-strength is located in 
discrete states where no ambiguity occurs, the s-strength 
is located in the continuum between 20 and ^50 MeV re- 
moval energy. In this region, also the correlated strength 
contributes, and affects the shape of the fitted WS mo- 
mentum distribution. The correlated strength has a mo- 
mentum distribution that falls much more slowly with 
increasing k than the Is QP strength; when fitting the 
sum with a WS parameterization, the resulting WS mo- 
mentum distribution extends somewhat too far in mo- 
mentum, i.e. it would have too small a radial extension 
in r-space. 

We have used the theoretical Is QP and correlated 
strength in the region used in Q for the determination 
of the Is momentum distribution to calculate a correc- 
tion to the WS parameterization fitted to the sum. We 
find that the extent of the wave function in k-space needs 
to be reduced by 11%. This modified WS-shape has been 
confirmed |32j by an independent analysis employing the 
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recent JLAB data of |5| and the correlated spectral func- 
tion of 0| . 

With these QP radial wave functions, and occupations 
renormalized to the one derived above from the high-g 
experiments, we can compute the QP density in r-spacc. 
As the radial wave functions fitted by Lapikas et al. re- 
fer to relative coordinate between the proton and the 
CM of the (A-l) system, we need to rescale the radial 
coordinate by a factor 11/12, with the corresponding ad- 
justment in height to conserve the normalization. 
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FIG. 2: Comparison of QP and correlated densities. 

The resulting QP density is shown in figs. 121- il The 
difference of the 12 C point density and this QP density 
yields the correlated density also shown in Figs. 12131 The 
'error bar' on the correlated density is not straightfor- 
ward to calculate due to the various adjustments that 
had to made in the analysis; ±20% for r < 3/m is prob- 
ably a realistic estimate. 



V. COMPARISON OF QP AND CORRELATED 
DENSITIES 

The first observation one can make from Figs. |2E1 con- 
cerns the fact that like in the theoretical analysis of Fig. Q] 
also the correlated density deduced from experimental 
data is significantly more concentrated towards the nu- 
clear interior than the QP density. One also observes 
that the correlated density in the nuclear interior gives a 
very significant contribution, of order 30%, of the central 
density, i.e. larger than one could have expected from 
the number of 20% or so of correlated nucleons. This 
explains why attempts to explain the total densities in 
terms of QP orbitals cannot be very successful. 

Figure [3] shows another important feature: The ex- 
perimental QP density at large r agrees perfectly with 




FIG. 3: Comparison of QP and correlated densities. The 
crosses indicate the tail of the density obtained when using 
the low-g QP occupation. 



the point density measured via clastic electron scatter- 
ing. Such an agreement should occur, as at large r — 
outside the range of the nuclear potential — the density 
is entirely given by the tail of the least-bound QP or- 
bit, the IP3/2 state in 12 C. More deeply bound states, 
or correlated nucleons with large removal energy, cannot 
contribute. 

The good agreement between QP density and point 
density at large r also confirms the correctness of our 
choice of QP occupations. Had we used the occupation 
derived from the low-g experiments and, as above, the 
shape of R(r) determined in Q from the fit to the world 
(e,e'p) data, we would have obtained the tail indicated 
in fig. |21 by the crosses. These are obviously significantly 
too low. This comparison thus provides an a posteriori 
justification of our procedure. 

In figs. l2!3l we also show the correlated density obtained 
by theory. Considering the above-mentioned uncertainty 
of the experimental result, we consider the agreement 
between theory and experiment as a good one. The size 
of the correlated contribution in the nuclear interior is 
very similar, the rapid fall-off of the correlated density 
at large r also agrees within our estimated uncertainty. 
The correlated density deduced from experimental data 
seems to contain more spectral strength in partial waves 
with I > than the theoretical one. 

The large contribution of the correlated density in the 
nuclear interior shows that the neglect of this correlated 
contribution in the standard IP calculations (e.g. all the 
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shell-model descriptions) is not justified. 

VI. CONCLUSIONS 

Starting from (e,e'p) data we have constructed the QP 
density for 12 C in coordinate space. The difference to 
the total density, obtained from elastic electron scatter- 
ing, provides the density distribution of the correlated 
nucleons. We find that it is significantly more concen- 
trated towards the nuclear interior. We also find good 
agreement with the theoretical calculation of the corre- 
lated density distribution. 

The large contribution of the density related to short- 
range N-N correlations, ^30% in the nuclear interior, 
together with the fact that the shape of the correlated 
density differs strongly from the QP density, explains 
the poor performance of QP wave functions in explaining 
many observables. Due to the shape difference, the short- 
coming of the neglect of the correlated contribution also 



cannot satisfactorily be 'compensated' by using effective 
quantities like effective charges, etc. 

As a side product, our analysis provides a solution to 
the puzzle raised in pj , the pronounced disagreement be- 
tween QP occupations derived from the low-q and high-g 
(e,e'p) experiments. We find that only the high-g occu- 
pation (suitably corrected for the correlated contribution 
not considered in ) is compatible with the independent 
information from elastic electron scattering. 
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